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Abstract 

Introduction. Functionally graded materials are of great use, because heterogeneity of properties enables to control the 
strength and rigidity of structures. This has caused great interest in the topic in the world scientific literature. The 
construction of solutions to such problems depends significantly on the type of boundary conditions. In this paper, we 
consider the equilibrium of a thin-walled circular cylinder whose mechanical properties change along the radius. 
Homogeneous boundary conditions were set on cylindrical surfaces that had not been considered before, the effect was 
on the ends. The mathematical formulation of the problem was carried out in the linear theory of elasticity in the 
framework of axisymmetric deformation. Expressions were constructed for the components of the stress-strain state of 
the cylinder, in which some coefficients were found from the solution to the resulting system of linear algebraic equations. 
Materials and Methods. The material of the cylinder was linearly elastic, the elastic modulus of which depended linearly 
on the radial coordinate. The basic research method was the asymptotic method, in which half the logarithm of the ratio 
of the outer and inner radii acted as a small parameter. Iterative processes were used to construct the characteristics of the 
stress-strain state of the cylinder. 

Results. Homogeneous solutions to the boundary value problem were obtained for a linearly elastic functionally gradient 
hollow thin-walled cylinder. An analysis of these solutions made it possible to reveal the nature of the stress-strain state 
in the cylinder wall. For this purpose, an asymptotic analysis of the solutions was carried out, relations for displacements 
and stresses were obtained. It was determined that those solutions corresponded to the boundary layer, while their first 
terms determined Saint-Venant edge effect similar to the plate theory. 

Discussion and Conclusion. The analytical solution to the equilibrium problem of a thin-walled cylinder inhomogeneous 
in radius constructed by asymptotic expansion can be used for numerical solution to a specific problem. For this, it is 
required to solve the obtained systems of linear algebraic equations and determine the corresponding coefficients. The 
resulting asymptotic representations provide analyzing the three-dimensional stress-strain state. The selection of the 
number of expansion terms makes it possible to calculate displacements and stresses with a given degree of accuracy. 


This analysis can be useful in assessing the adequacy of applied calculation methods used in engineering practice. 
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AHHOTalna 

Beedenue, DyHKMoHaIbHO-rpamqMeHTHble MaTepHasIbl HAXOLAT OOMbIIOe IPHMeHeHHe, T.K. HEOMHOPOAHOCTh CBOMCTB 
HO3BOACT YIPaBIATh MPOUHOCThIO MW %KCCTKOCTHIO KOHTPYKIMH. ITHM BbI3BAH OOIbIIOM HHTepec K WaHHOM Teme B 
MUpoOBOM Hay4dHoH WuTepatype. locrpoeHue pellieHua TakuX 3ajad CyNIeCTBCHHO 3aBHCHT OT THIa TpaHH4HBIX yCJIOBHI. 
B Hactosileti padoTe paccmaTpHBaeTca paBHOBeCHe TOHKOCTeHHOrO KpyrOBOrO WWIMHApa, MexaHWueckHe CBOMCTBa 
KOTOPOTO 34MCHAIOTCA BAONb paguyca. Ha WuMHApH4yecKkHXx MOBePXHOCTAX 3aaHbI OMHOPOAHbIe TPaHWGHble YCOBUA, 
KOTOpbIe WO ITOTO He paccMaTpHBAaJINCh, BO3elCTBHe OKa3bIBaeTCA Ha Topulax. MarematTw4eckaad NocTaHoBKa 3aa4qu 
OcyllecTBIAeTCA B JMHeMHOM TeopHu yipyrocTH B paMKax OCeCHMMeTpHYHON Aedopmaiuuu. B padote mocTpoeHsI 
BbIP@KEHHA JI KOMIOHCHT HalipsKeHHO-eOpMMpoBaHHOTO COCTOAHUA IWIMHApa, B KOTOPBIX HeKOTOpbIe 
KOSPPUIUMCHTHI HAXOAATCA M3 PeMICHHA MOWYICHHOM CHCTeMBI JIMHeMHBIX alreOpav4eckux ypaBHeHHii. 

Mamepuaaoi u memooot. Martepual WunuHypa ABIIAeTCA JMHeEMHO yIPyruM, MOLyJIb ypyrocTu KoTOporo JIMHeiHHO 
3aBHCHT OT payvaIbHOM KOOpAMHAaTbI. OCHOBHBIM MeTOJOM HCCIeOBaHHA ABJIAeCTCA ACCHMITOTH4eCKHH MeTOL, B 
KOTOPOM B KayecTBe MaJIOrO MapaMeTpa BbICTyHaeT MOMOBHHA JOrapHdMa OTHOINeCHHA BHeINHerO UH BHYTpeHHero 
panuycos. Jia NocTpoeHua XapakTepHcTHK Hallps2%KeHHO-eOpMMpoBaHHOrO COCTOAHUA IMIMHApa TMpHMeHeHbI 
UTepallMOHHBle IIpOleccsl. 

Pezynomamet ucciedosanua. Jina sMHewHO-yupyroro (dyHKIMOHaIbHO-rpagqWeHTHOTO MONOrO TOHKOCTeHHOTO 
IWIN Apa WOUYYeCHEI OTHOPOAHbIe pellieHHa KpaeBoli 3ayaun. AHasIH3 3TUX PeIICHuM MO3BOIMeT PacKPbITh XapakTep 
Hallps.KeHHO-edOpMHpOBaHHOTO COCTOAHHA B CTeHKe IMIMHApa. C ITO Lebo MpoBeseH ACHMITOTHYeCKHH aHasH3 
pelleHuii, MOJyUeHbI COOTHOMICHHA JIA MepeMelleHHi MU HallpsyKeHHM. YcTaHOBJICHO, UTO ITH PeIIICHHA COOTBETCTBYIOT 
MOrpaHHYHOMy CJIOIO, Ip ITOM HX MepBble YWIeHbI OMpeeuAOT KpaeBont 9dext Cen-Benana, aHasOrM4Hblli Teopuu 
TMT. 

O6écyotcoenue u 3axu04enue. TloctpoeHHoe Cc MOMOLIbIO ACHMIITOTHYECKOLO pa3zIOXKeHHA aAHAIMTMYeCKOe pelleHue 
3alauH O paBHOBeCHH HeOHOpOAHOrO 0 paswycy TOHKOCTeHHOrO IMIMHEpa MO#KeT ObITb HCIONb30BaHO IIA 
YMCJICHHOrO pelieHHA KOHKpeTHOM 3aqauu. Jina 9TOrO HYyKHO pelIHTb TMOyYeHHbIe CHCTeMBI JIMHeMHBIX 
anreOpawuecKux ypaBHeHHii WM OlpeseuMTb coorBeTcTByroulMe Kosddunnents. TlomyaeHHple acuMmToTH4eckHe 
IIpeJicTaBICHHA TO3BOUAIOT aHaIM3HPOBaTb TPeXMepHOe HalipsxKeHHO-ed@opMupoBaHHoe cocTosHHe. Bsroop 
KOJIM4YECTBA YICHOB pa3zIO2XKCHUA MO3BOIAeCT PACCUHTATb WepeMelCHHA H HalIpMKCHHA C 3a/JaHHOM CTeMeHbIO TOUHOCTH. 
OTOT aHaiv3 MOxKeT OBITh MONe3eH MPH OWeHKe aleKBAaTHOCTH UpPWKaqHbIX MeTOOB pacueTa, IPHMeCHAeMBIX B 


MHKeHeEpHOH lpakTuKe. 


KsoyeBbie CioBa: JMHelHadt Teopua yupyrocTu, PyHKUHOHANIbHO-rpaqMeHTHbIM MaTepHasl, TOHKOCTeHHbIM MOJIbIM 


UWHJIMVHAp, OMHOPOAHbIe pelwleHHusA, nMorpaHvuHEit con, BapHaliMOHHbIit TIPHHUAIT 


BaarogxapHocrnu: aBTOPbI BbIpaxKaroT OaroapHOcTh pedakWHu U peleH3enHTy 3a BHHMAaTeJIBHOe OTHOLMICHHE K CTaTbe UH 


TIPCAIO7KCHHA, KOTOPbIe MO3BOJIMJIM HOBbICHTb Ce KavyecTBO. 
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Introduction. Functionally graded materials are widely used in various constructions. Due to the dependence of 
mechanical properties on coordinates, it is possible to control the stress-strain state (SSS) of parts. An example of using 
such inhomogeneity is a cylinder whose mechanical properties depend on the radius. In this case, the cylinder may be of 
interest as a separate structure, or as being a division subring of a compound body, e.g., connecting two media with widely 
different properties. When calculating the SSS of a thin-walled cylinder, some applied theories can be used, an assessment 
of their adequacy. Especially in the case of inhomogeneous properties, it can be carried out using computer modeling or 
asymptotic analysis based on a three-dimensional formulation. The latter determines the relevance of this study. 

A number of studies have been devoted to the investigation of the SSS of hollow cylindrical bodies within the 
framework of the linear theory of elasticity. In [1, 2], the mechanical behavior of a radially inhomogeneous cylinder in a 
three-dimensional formulation was studied on the basis of the spline collocation method and the finite element method. 
In [3], the SSS of a cylinder whose properties depended on the radius loaded with uniform internal pressure was described. 
In [4], an analytical study was carried out for a functionally graded piezoelectric cylinder. In [5], an exact solution was 
constructed to a radially inhomogeneous hollow cylinder with exponential Young's modulus, with constant Poisson ratio 
and power Young's modulus. In [6, 7], an analytical solution to the axisymmetric thermoelasticity problem for a 
continuous cylinder was obtained using the direct integration method when the coefficient of linear thermal expansion 
was an arbitrary function of the radius. In [8], a general asymptotic theory of a transversally isotropic homogeneous 
hollow cylinder was developed. New groups of solutions were obtained for a transversally isotropic homogeneous 
cylinder. The comparison of the constructed solutions to the solutions constructed using applied calculation methods was 
given. In [9, 10], some boundary value problems of elasticity theory were studied for a functionally gradient isotropic and 
transversally isotropic (the plane of isotropy was perpendicular to the axis) cylinder, in the case when the elastic modules 
were arbitrary continuous functions of the radius of the cylinder. In [11], an analysis of the bending deformation problem 
for a radially inhomogeneous cylinder was carried out. The analysis of the above papers shows that not all types of 
boundary conditions on cylindrical surfaces have asymptotic representations of solutions. 

In this article, on the basis of an asymptotic analysis of three-dimensional equations of elasticity theory, the features 
of the SSS of a thin-walled cylinder whose properties vary linearly along the radius were studied. In this case, the inner 
border was fixed in the axial direction and was free in the radial direction. 

This goal was achieved using several steps: asymptotic integration of differential equations and the construction of 
homogeneous solutions; derivation of formulas for the components of the displacement vector and stress tensor; 
consideration of boundary conditions on the face surfaces. 


Materials and Methods. Radially inhomogeneous hollow thin-walled cylinder 
T= {r E iG ; Fal be [0; 2r], ze [=i Ly } is considered in a cylindrical coordinate system with the origin on its axis. The 
problem of its equilibrium, in the case of fixing cylindrical surfaces along the axis and zero normal stresses, is solved in 


an axisymmetric formulation under the action of stresses at its faces. 


The boundary value problem consists of the equilibrium equations [8]: 


Oo 0o 19 oe O45 


rT p 'Z 4 — 0) ‘ 1 

or Oz r ) 
09, i 09. re oO, 0, (2) 
or Oz r 


where 6,,,0,,. O4)>O,, — components of the stress tensor. 


rr? 


Defining relations [8]: 


s, = (2642) vi +) (3) 
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0, =(26+2) aH vi rte) (4) 
7 Oz or r 
u Ou. Ou 
6,, =(2G+A)— +A “4+ |, a) 
” ( ) r (= Oz a 
S. -6{& aa (6) 
. oz Or 


Here, uv, =u,(7,z), u, =u,(r,Z) — components of the displacement vector. 
Lamé parameters vary linearly along the radius: 
G(r) =G.,r, ArH dr, (7) 
where G,,A, — constants. 


After substituting (3)-(7) into equations (1), (2), the dimensionless system of equations takes the form: 


O2u Ou Ou, Ou, O2u 
(2G, +,)| —? +e— |+e(G, +A, Jew +1,62e? —+G,62e% —* —2G,6u, =0, (8) 
op? Op Opes 0§ 08? 
O2u Ou Oru Ou (ol 
G,| —= +e— |+(2G, +A, )e?| e% —++e" — |+(G, +A, )ee? —*2 =0. (9) 
op? Op 0S? 0§ Opog 
Here: 
1 r Zz ; ‘ : 1 Lh ; ‘ 
p=—lIn a new dimensionless coordinates; ¢ = a — | —in the case of wall thinness, small parameter; 
é Nt % q 
ii, Uu,. u. el G.r, : 
m=Vnn, pe[-bl], 6e[-Gl], /=2;u,=+, w=, Y= , G, =—*; G, —some parameter having the 
% % % G, G, 


dimension of stress. 


Let us consider a problem in which homogeneous boundary conditions are set on the lateral surfaces of the cylinder: 


Melo =O (10) 
Spl, =0- (1) 
Stresses are applied to the faces of the cylinder: 
Orefe uy = 4s (P)> (12) 
Oes|. 4, =Ls(P)» (13) 
(St). 
oO or O,. = Pr: O,. = ox dimensionless stresses 
pp G, 7 PE G, 2s . 


Stress vector components 1,,(p),t,,(P), (s=1;2) satisfy the equilibrium conditions. 


To construct homogeneous solutions, we seek the components of the displacement vector in the form: 


u,(psE)=u(p)ews, u,(p:5)=w(p)er. (14) 

Substituting representations (14) into the system (8)—(11), we obtain: 
(2G, +A, )(u"(p) +eu' (p)) + eae ((G, +A, )w'(p) +A ew(p)) + €2G, (a2e2 —2)u(p) =0, (15) 
G,(w"(p)+ew' (p))+(2G, +A, )e? (aeu(p)+ a2e2w(p))+e(G, +A, )aeru' (p) =0, (16) 
wo, = 9, (17) 
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=0. 


p=+1 


[ (2G, +o )u'(p)+ Ero (u (p)+ aew(p)) 


(18) 


We investigate boundary value problems (15)-(18) at se +0. To solve (15)-(18) at ¢ +0, we use the asymptotic 


method [9-13]. 


Nonzero solutions (15)—(18) correspond to the third iterative process, the components of the displacement vector are 


searched for in the form of expansions over a small parameter: 


i() Selec (peace) a), 
w®) (p) = (i (p) + ew, (p)+---), 
a=e!(B, +eB, +---). 


After substituting expansions (19) into equations (15)—(18) for terms of the first order, we have: 
(2G, a5 Ay U5 (P) +By (G, +o) Way (P) +G, aU39 (p) = 0, 
GW (p) +B, (G, + Do) Miho (p) + (2G, +o )B3 Wo (p) =0, 


=0, 


p=tl 


((2G, +A )usy (p) +A BoW5o (p)) 


Following [13], the spectral problem (20)—(23) corresponds to a potential solution for the plate. 


Wo (p) 


p=+1 . 


Thus, the solutions are presented in the form: 


a) UG) (p38) = exT, (-2G,B3, sinBy, sn(up)*+0(@))exp{ *(B +eBi, +e) s| . 


un (p38) - EDT, (2G, dx Si Bo, £05(Bp) +O(6))<exp{ * (Bu +B, + es) 6] : 


Here, B,, is the solution to the equation: 
cos Bo, =0. 


The stresses corresponding to solutions (24), (25) have the form: 


a” - DT, (—4G2B3, sin Box £05(Bup) +0(6))<exp{ * (Bu +B, +++) s] ’ 


os” = 3 T, (4G38;, sin By, sin(yp)+0(6))ex0| “(Bo +e€B,, + oo) | > 
kal € 


= 


3:1 < 300 1 
Gy y= 2 [, (4G28;, sin Bo, £05 (BP) +O(@)}exo{ “(Pa +eB,, +++) s| 
oi” - (6). 


b) uO) (p; é) - ey F, (2G,B3, cos By, cos (Bo:P) ss O(e)) “ exp( “(Pu +é€B,, ++ JB) 


i=l 


ue (p; 5) = ey F, (2G,B3, COS By, sin(Bo;P) + O(e)) “ exp( 4 (Bo, + eB, ++ ) s| , 


Here, B,, is the solution to the equation: 
sin B,; = 0. 


The stresses corresponding to solutions (31), (32) have the form: 


O82) = LF (—4G38}, cos By, sin (By,P) +O(c))es9{ “(By +éB,, +36} 


of = LF, (46283, cos B,, c05(Bup)+0(€))exP( 2 +eB,, +6} 


(19) 


(20) 
(21) 


(22) 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


(35) 
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0° = 5 F (4638, cosB,, sin (Bp) +0(6))err( (Bs —— 8} ” 


of” =O(e). (37) 
The general solution (15)—(18) will be a superposition of solutions (24), (25), (31), (32): 


u, (p:&) =) 7, (-2G,Bi, sin Bo, sin (By,p) + O(e))x 
xexp( M(B +eB, + ++) g}robF (2G,B3, cos By; cos(By,P) + 0(@))exp( (Pn +éB,,+ oo) é} (38) 


Ue (p;6) = ey T, (2G,B%, sin, cos(By,P) + O(e))x 


> 
un 


«exp( 1B +B, +e) eS (2G,B3, cos By, sn(yp)+0(@) exp 2(By, +éB,, +98), (39) 


i=l 
Solutions (24), (25), (31), (32) have the character of a boundary layer. When moving away from the faces, solutions 


(24), (25), (31), (32) decrease exponentially. 
To determine constants 7,, F,, we use Lagrange variational principle. The variational principle takes the form [8]: 


i 


oy AL 
2 (6. fis OM, * (ox hs Jou, | Loe =0. (40) 
Substituting (24-36) into (40), we have: 
EM Tug = Paz: (41) 
$0,Fo = pl a 


Here: 


M x = 16G5B5; ok (i — Box i‘ sin Bo, sin Bo; sin (By, “Pao Bafa g{ Oa =Pul} (npu j# k) 
, € 
M ,, =16G;B3, sin? By, [exe ei 24) (npu j =k) 
1 l 
Po; = 2G}; sin Bo, it. (p)cos(B,,p)-t, (0)sin(B,p))do-evo{- a 
1 
Q;, = 16G5B5 Bo, (B,, — Boi y" COS By ; COS Bo, sin (By, —Bo: \e*[- Gest), Saal (npu i# j) 


2p, 2B, / 
Q, = 166585; cos’ By; feo{ 2). e( 74) (npu i=j ) 
, €& €& 


| oj! 
Poy) = 2G,B3; COSBy, ue (p)cos(B,,p) +12, (p)sin(B,,p))dp- cxn(- re }: 


+5 (tm (p)cos(B, ,P) —t,,(p)sin (B,,p))dp cx © 


Jj 
€ 


| . By jf 
+] fie (p)cos(B,,P) +t, (p)sin(B,,p))dp exe{"“}} 
TTF ely #4; 

F, = Fy +éF, Hees 
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Constants T,,,F,, (p =1, 2,...) are found from systems of linear algebraic equations (41), (42), analogous to which are 
studied in [13]. 

Research Results. Tp the article, in an axisymmetric formulation, the solution to the problem of linear elasticity theory 
for a functionally graded hollow thin-walled cylinder, whose properties vary in thickness according to a linear law, was 
considered. Homogeneous cross boundary conditions were set on the lateral surfaces of the cylinder, and a stress vector 
was set at the faces. The constructed homogeneous solutions satisfied boundary conditions on cylindrical surfaces. For 
their construction, an asymptotic approach based on the expansion by a small parameter characterizing the relative 
thickness of the cylinder was used. To account for inhomogeneous boundary conditions at the faces, systems of linear 
algebraic equations similar to those studied in the literature were obtained. It was shown that the constructed SSS solutions 
had a boundary-layer character, which corresponded to the edge effect similar to the theory of inhomogeneous plates, 
which bears the name of Saint-Venant. 

Discussion and Conclusion. Usually, when studying the SSS of thin-walled structures, applied calculation methods 
are built that reduce the dimension of the problem. In this regard, the task of determining the range of geometric and 
mechanical parameters in which these methods give acceptable accuracy is critical. The solutions to three-dimensional 
equations constructed in the work on the basis of asymptotic analysis make it possible to assess the adequacy of such 
applied theories with a predetermined accuracy threshold. In addition, these solutions can find application in the 


evaluation of numerical solutions to problems for structures with functionally graded materials. 
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Kouqbauxm unmepecose: aBtTop 3a0BiAeT OO OTCYTCTBHH KOH(JIMKTa HHTepecos. 


Aemop npouwumas u OdObpul OKOHYAMeEbHbI BapUaHM pyKoNucU. 


